CTEQ/408 
MSUHEP-50319 
April, 1995 



Effects of QCD Resummation on Distributions of 
Leptons from the Decay of Electroweak Vector Bosons 



Csaba Balazs M, Jianwei Qiu ^ and C.-P. Yuan ( a ) 

(a) Department of Physics and Astronomy, Michigan State University 
East Lansing, MI 48824, U.S.A. 
(b) Department of Physics and Astronomy, Iowa State University 
Ames, IA 50011, U.S.A. 

Abstract 

We study the distributions of leptons from the decay of electroweak vector 
bosons produced in hadron collisions. The effects of the initial state multi- 
ple soft-gluon emission, using the Collins-Soper resummation formalism, are 
included. The resummed results are compared with the next-to-leading-order 
results for the distributions of the transverse momentum, rapidity asymmetry, 
and azimuthal angle of the decay leptons. 



1 Introduction 



With the discovery of the top quark , the electroweak symmetry breaking mech- 
anism remains one of the major mysteries of particle physics today. Unfortunately, 
the precision low energy data have told us very little about the scalar sector, i.e. the 
electroweak symmetry breaking sector, of the Standard Model (SM). From CERN 
LEP data we learned that the mass m# of the Higgs boson has to be larger than 
about 60 GeV [[J . However, the precision low energy data do not exclude the possi- 
bility for m,u to be in the order of 1 TeV. At the moment, one of the largest theoretical 
errors in analyzing radiative corrections to low energy data comes from the predic- 
tion of the fine structure constant a evaluated at the Z-boson mass scale, due to the 
less precise low energy e~e + — *> hadrons data 0. Without further improvement in 
the determination of a(Mz), a more precise way to test the SM is to have a better 
measurement of Myy. If M^y is measured within 40MeV and the mass of the top 
quark within 4 GeV, then m# can be constrained within a couple of hundred GeV 
!|. To reach such an accuracy in the measurement of My/ at hadron colliders, we 
have to know the kinematics of the W^-boson well. Since W + decays into a charged 
lepton £ + and a neutrino U£, the kinematics of the W + cannot be accurately known 
because of the missing momentum carried by vg_. It is therefore desirable to have a 
good prediction on the kinematics of £ from the decay of W ± . 

The well established fact that the transverse momentum Qt distribution of the 
IV ± -boson cannot be described by the next-to-leading-order (NLO) perturbative cal- 
culation in low Qt region |J implies that the transverse momentum p T of the lepton 
^ cannot be accurately predicted by the NLO calculation, especially for p T ~ M w /2 
(mostly with low Qt) where the data dominate. We must resum the effects of the ini- 
tial state multiple soft-gluon emission to predict the distributions of the leptons from 
the decay of the vector boson V (= 7, W or Z) produced in hadron collisions!] In 
this paper, we adopt the Collins-Soper formalism ||, and closely follow the notation 

1 The analytic results presented in the paper also apply to the non-standard weak gauge boson(s) 
(such as Z') present in any extended gauge theory. 
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used in Ref. [[F] to resum the multiple soft-gluon effects to the transverse momentum, 
rapidity asymmetry, and azimuthal angle distributions of the decay product leptons. 



2 The Resummation Formalism 

To obtain the resummed results, we use dimensional regularization to regulate 
the infrared (IR) divergencies, and adopt the canonical-75 prescription to calculate 
the anti-symmetric part of the matrix element in n-dimensional space-time. Q The 
infrared- anomalous contribution arising from using the canonical-75 prescription was 



carefully handled by applying the procedures outlined in Ref. |L0|] for calculating both 
the virtual and the real diagrams.^ 

The kinematics of the vector boson V (real or virtual) can be expressed in the 
terms of its mass Q, rapidity y, transverse momentum Qt, and azimuthal angle <pvi 
measured in the laboratory frame. The kinematics of the leptons from the decay of 
the vector boson can be described by the polar angle 9 and the azimuthal angle 0, 
defined in the Collins- Soper frame [|ll|], which is a special rest frame of the l/-boson 



12| . The four-momentum of the decay product fermion in the lab frame i 

ru Q /V N 



where 



+ sin cos + sin #3^01^ + cos # , (1) 



g M = (Mt cosh y, Qt cos 4>, Qt sin 0, Mt sinh y) , 

X ^-^M- T ( q+n " + q ~ n "-f q ^ 
Z " = W T {q+n>1 ~ ^ ' 

Y » = e^ a ^Z a Xp. (2) 



Here, q± = -±(q°±q 3 ), M T = ^Q 2 + Q 2 T , n v = ^(1, 0, 0, 1) and n v = ^(1, 0, 0, -1 



2 In this prescription, 75 anticommutes with other 7's in the first four dimensions and commutes 
in others [§, i_ 

3 In Ref. Jl0| | the authors calculated the anti-symmetric structure function F3 for deep-inelastic 
scattering. 

4 Our convention is that — (q° . q 1 , q 2 , q 3 ), Q = y/q 2 , and Qt = \l {q 1 ) 2 + (<Z 2 ) 2 - The total 
anti-symmetric tensor e 0123 = — 1. The proton beam direction is assigned to be the positive z-axis. 
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To obtain the fully differential cross section of the vector boson production and 
decay for all values of Qt, we need the resummation formula [|7|: 

/ da(AB -> U')X) \ 1 Q 2 



dQ 2 dy dQj, dcf)y d cos 9 1 



9Qir 2 S (Q 2 -M 2 ) 2 + M 2 T\ 



I. 



d 2 b W jk (K, Q, x A , x B , 9, 0) F^ib, Q, x A , x B ) 



\(2tt) 2 
+ Y(Q T ,Q,x A ,x B ,9,(f))\ 



(3) 



Here Wj k is 

W jk {b, Q, x A , x B , 9, 4>) = exp {-S(b, Q)} \ V jk \ 2 
x { [{C ja <g> f a/A ) (x A ) (C- kb ® f b/B ) (x B ) + (C- ka <g> f a/A ) (x A ) (C jb ® f b/B ) (x B ) 
x(g 2 L + gR)(fl + fn)(l + cos 2 9) 

+ [(Cja ® fa/A) (X A ) (C- kb ® / b /s) (x B ) ~ (C ka ® / a/A ) (^a) (C^fe <g> / 6 / S ) (^fi) 

x G£ - fl!) (/£-/£) (2 cos 0)}, (4) 
where ® denotes the convolution and is defined by 



(C ja <8> /a/^) Oa) 



1 <&4 



. /a/A (U,t*) C ja \-z-,b,fJ, 
*A ?A VsA y 



(5) 



and the coefficients are given by 



J Cabibbo — Kobayashi — Maskawa matrix elements for V = W ± 
jk= \6 lk forV = Z°, 7 



(6) 



In the above expressions j represents quark flavors and k stands for anti-quark flavors. 
The dummy indices a and b are meant to sum over quarks and anti-quarks or gluons. 
Summation on these double indices is implied. 
The Sudakov form factor S(b, Q) is defined as 



5(6, g) 



rQ 2 dp 2 


In 


(Q 2 \ 


lbl/b* fL 2 







(7) 



The A, B functions and the Wilson coefficients Cj a , etc., were given in Ref. [0]. After 
fixing the renormalization constants C\ = bo = 2e _7B and C% = 1, one can obtain 
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A®, A® and B^ from the Eqs. (3.19) to (3.22) of Ref. §.Q ( 7e is the Euler 
constant.) After choosing /z such that fib = C3 = 2e _7E , the Wilson coefficients 
for the parity-conserving part of the resummed result are greatly simplified from the 
Eqs. (3.23) to (3.26) of Ref. as 



C$ = 5 ]k !^(l-z)+ 1 -(n*-8)5(l-z)} and = -<!-:). 



(8) 



Following the procedures given in Ref. |L0[ for handling the 75's in n-dimensional 
space-time, we find that the same Wilson coefficients also apply to the parity- 
violating part of the resummed result. 

In Eq. (|^), the impact parameter b is to be integrated from to 00. However, for 
b > b max , which corresponds to an energy scale less than l/b max , the QCD coupling 
a s becomes so large that a perturbative calculation is no longer reliable. f] Hence, 
the non-perturbative function F NP is needed in the formalism, and generally has the 
structure 



Fjl P (b,Q,Qo,XA,x B ) = exp -In ^^^J h x (b) - h j/A (x A ,b) - h~ k/B 



[x B ,b) 



(9) 



where hi, hj/ A and %/ B cannot be calculated using perturbation theory, so they must 
be measured experimentally. Furthermore, W is evaluated at 6*, with 

b 



6* = —,= 

yfl + (6/6, 

such that K never exceeds b max - 

The F-term in Eq. (|3[) is defined as 

= f 1 dU 

Jxa £,A J kb £,B n= i 



(10) 



Y(Q T ,Q,x A ,x B ,9, 



1 d^ B 



E 



Ots{Q) 



71 



N 



Xfa/A^A] Q) RdP(QT, Q, ZA, ZB, 0, 4>, Q) fb/B^B] Q), 



(11) 



in which the functions only contain contributions which are less singular than 
Qt 2 x (logs or 1) as Qt — * 0. Their explicit expressions for pp — > V(— >• ll')X are 
given in the Appendix. 



J For instance, we obtain = 4/3 and B^> = -2. In our numerical results we also include 
and flW. 



We use b„ 



0.5 GeV in our calculation. 



3 Numerical Results 



In this paper, we only give our numerical results for pp — > W + {— > £ + vi)X at the 
Fermilab Tevatron with = 1.8 TeV. The CTEQ3M parton distribution functions 
(PDF's) are used along with the non-perturbative function [13] 

F NP (b, Q, Q Q , x A , x B ) = exp | - 9l b 2 - g 2 b 2 In (J^J - 9l g 3 b In (100x A x B ) j , (12) 

where g x = 0.11 GeV 2 , g 2 = 0.58 GeV 2 , g 3 = -1.5 GeV" 1 and Q = 1.6GeV.Q To 
consistently compare the distributions of the leptons in NLO and resummed calcula- 
tions, we have used exactly the same PDF's, QCD and electroweak parameters, etc., 
for calculating the NLO results.^ Furthermore, we have applied the kinematic cuts 
p T > 25 GeV, $t > 25 GeV, and Qt < 20 GeV. These cuts are similar to those ap- 
plied by the CDF group in the measurement of the asymmetry in the lepton rapidity 
distribution from H^-boson decays.0 

The transverse momentum distributions of the charged lepton p T are shown in 
Fig. [I] for NLO and resummed calculations. We note that in these results a Breit- 
Wigner resonant width has been included, cf. Eq. (|3|). In the vicinity of p l T = Mwjl 
("Jacobian peak") the NLO calculation is ill-defined because its amplitude blows 
up as Qt — > 0. The effect of the initial state multiple soft-gluon emission on the 
distribution of p T is to widen and smoothen the "Jacobian peak" and therefore make 
it more challenging to accurately extract My/ from the p T distribution. 

Recently, the Fermilab CDF group measured the asymmetry A y t in the rapidity y e 
distribution of the charged lepton (£) from the decay of W + — > l^vg (or W~ — > i~vt) 



[]nj| and proved this measurement to be particularly sensitive to the slope of the ratio 

of u- to d-quark parton densities inside the proton [H5, [T7J. Define the asymmetry 

in the lepton rapidity distribution as 

da/dyH/ > 0) - daldy\y l < 0) 
yl da/dy i (y?>0) + da/dyt(y e <0)' 1 ' 



7 These values were fit for CTEQ2M PDF, and in principle should be refit for CTEQ3M PDF. 

8 Our NLO results agree with those in Ref . Hj . 

9 The requirement of Qt < 20 GeV in our calculation is approximately equivalent to cutting out 
the events in which the transverse momentum of the net hadronic activities reconstructed from the 
calorimeter cells within the pseudo-rapidity range of ±3.5 is larger than 20 GeV. 
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27.5 30 32.5 " 35 37.5 40 ~ 42.5 ~ 45 

P £ T (GeV) 



Figure 1: Transverse momentum distribution of the charged lepton p T for NLO 
(dashed) and resummed (solid) calculations. Resumming the initial state multiple 
soft-gluon emission has the typical effect of smoothening and widening the Jacobian 
peak (at = M w /2). The NLO distribution is singular and ill-defined near M w /2. 




Figure 2: The asymmetry A y t in the lepton rapidity distribution as a function of 
y for NLO (dashed) and resummed (solid) calculations. They differ the most in the 
large rapidity region (\ye+\ > 1). The experimental data were obtained from Ref. [15] 
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Figure 3: The distribution of the difference in the lepton azimuthal angles near 
the region A0 ~ tt. The NLO (dashed) distribution is ill-defined at A0 = it and 
is arbitrary around it. The resummed (solid) distribution gives the correct angular 
correlation of the lepton pair. The distribution has a similar peak for A<p ~ — ir. 

which is commonly known as the lepton charge asymmetry. In Fig. 0, we show A y i 
as a function of y l for NLO and resummed calculations. As indicated in the figure, 
they differ the most in the large rapidity region (\ye+\ > 1). Recall that as Qt — > 
the NLO Qt distribution becomes singular, but the resummed result remains finite. 
The rapidity y distribution of the jy-boson in the NLO calculation is not singular, 
and we expect that after integrating out the complete phase space for Qt (that 
is, without imposing any kinematic cuts) the NLO and the resummed calculations 
should predict the same y distributions. We have explicitly checked that this indeed 
is the case. However, in Fig. |2] some kinematic cuts (as described at the beginning 
of this section) have been applied to our calculations. In the rest (Collins-Soper) 
frame of the H^-boson, the decay kinematics of the lepton £ is identical for both the 
NLO and the resummed calculations because the decay leptons do not involve strong 
interactions. Since the iy-bosons have different kinematic distributions (e.g., Qt 
distributions) in these two calculations, the resulting lepton kinematic distributions 
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(e.g., y distributions) in the laboratory frame are different. The two y e distributions 
differ the most in the large y region, where the typical y is large, since the effects 
of soft-gluon emission become more important, close to the boundary of the phase 
space. 

Another interesting observable to test the QCD theory beyond the fixed-order 
perturbative calculation is the measurement of the difference in the azimuthal angles 
of / and I' from the decay of V. In practice, this is better measured for Z — > £ + £r . For 
the sake of argument, we show in Fig. |3] the difference (A0) in the azimuthal angles 
of £ + and U£ measured in the laboratory frame for W + — > £ + v^ and calculated in NLO 
and resummed approaches. As clearly indicated, the NLO result is ill-defined in the 
vicinity of A<p ~ ±7r, where the multiple soft-gluon radiations have to be resummed 
to obtain physical predictions. Therefore, the transverse mass distributions of £ + -f£ 
pair for NLO and resummed calculations are also different ,0 and in principle only the 
resummed results can sensibly predict the distributions of the leptons for a precise 
measurement of My/. 

In conclusion, we found that the distributions (pj>, y and A0) of leptons are 
different in NLO and resummed calculations. For a better measurement of My/ and 
A y i , the effects of the initial state multiple soft-gluon emission have to be considered 
in hadron collisions. The more detailed phenomenological studies will be presented 
elsewhere. 
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The transverse mass of £ + -vi pair is defined as rojf 




9 



Appendix 



Let us define the qq'V and the IV V vertices, respectively, as 



Hn - 75) + 9r(X + 75)] and i-f^ [f L (l - 75) + f R (l + 75)] 



(14) 



For example, for V = W + ,q = u, q' = d, I = u e , and V = e + , the couplings 
gl = fl = G F M^/y/2 and g\ = j\ = 0. (G F is the Fermi constant.) In Eq. dTTj) , 
for N = h 



R 



(i) 

(lb 



16 I 



nQ 2 



(91 + g R m + r R )Rf + (Hi - 9km - f R )K , (15) 



where the coefficient functions Rf b are given as follows: 



f = r fk C + 7 ^±^1 5( s + t + u- Q 2 ) 

s 

f = r i~ k A 3 + K+ ^ u) 6(s + t + u- Q 2 ] 



Q 



A0 + A2 + -£-Ai 

VT 



M 2 ' 



2Q 2 K-(t,u] 



A 4 



Rf = r 9j C 







Q 2 T \M T Q T M T U + {t,u) 4 j' 

Q 2 Q\ H+(u,s) 



5(s + t + u-Q 2 



< W+(", -») [A + M + g_(<?-«)' + tt.(M)^ 



7Z + (u, s 



Qi 



7Z + (u, s 



X 



Q_ 



5(s + t + u~Q 2 ) 
Q (2u{Q 2 -s) 



u s 

2 



M T \ K + (u,s) 



1-1 



A, 



2Q 2 



Q T M 7 



2s{Q 2 



K+{u,s) 



+ 1 



Aa 



with 



r^+M 6{s + t+u _ Q 2 ) _ 2 6{1 _ Za) 5(1 _ Zb) 

LJrp S 



In ( ^ 

\Qtj 



(16) 



5(1 - *a) 



'1 + z 2 



1 - 2B, 



5(1 " Zfl) 



'1+4' 

1 - . 



3 
2 

(17) 



and 



r9j = £ 5(a + 1 + u _ Q2) _ r 2 + (1 _ ^ )2 j s{1 _ Zb) \ (lg) 



Q 2 



T 
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where lZ±(t, u) = (Q 2 — t) 2 ± (Q 2 — u) 2 . The angular dependence is described by the 
functions 

£ = 1 + cos 2 9, „4o = 7^(1 — 3 cos 2 9), Ai — sin 29 cos 0, A2 = ^ sin 2 9 cos 20, 
„4 3 = 2cos#, = sin 9 cos 0, (19) 

of which ^4 3 and At are odd under parity operation. 
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